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Drift  Ballooning  Modes  in  Multlpoles 

MARSHALL  N.  ROSENBLUTH  and  B.  COPPI 

University  of  California,  San  Diego 
and  General  Atomic  Division/General  Dynamics  Corporation 


It  is  well  known  that  In  the  presence  of  unfavorable 
curvature  the  drift  mode  is  unstable  even  for  zero  Larmor 
radius.   We  consider  an  electrostatic  perturbation  of  a 
multlpole  equilibrium  with  variable  curvature  in  the  limit 
that  the  frequency  is  high  compared  to  the  frequency  of  ion 
longitudinal  motion  but  low  compared  to  the  electron  fre- 
quencies.  In  this  limit  the  electrons  adjust  to  a  Boltzmann 
distribution  along  the  field  lines  while  the  ions  move  with 

■p  y  -p 

the  — p—  drift.   In  addition  an  electron  Landau  damping 

B 
will    exist  for  slow  particles  satisfying  oo  -   na)„  where 

CD„  is  the  frequency  of  the  longitudinal  motion.   This  will 

be  destabilizing  if  oo  <  ka.  —  ^  v.,  .  T  /T.  .   For  modes 

localized  in  the  neighborhood  of  the  unstable  curvature  this 

will  be  true.   We  find  that  the  condition  for  the  existence 

of  such  localized  ballooning  instabilities  is  given  by 


fkV)  /  1  dB  1  dn  ^ 

^K  a^j  ji     B  (jr  n  dr  T^  >  -^ 


where  £,   is  the  length  of  the  region  of  unfavorable  curvature, 
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Low  Frequency  Limit  of  Interchange  Stability 

MARSHALL  N.  ROSENBLUTH 

University  of  California,  San  Diego 
and  General  Atomic  Division/ General  Dynamics  Corporation 


For  simplicity  we  consider  an  axisymmetrlc  torus,  e.g., 
a  multipole  with  coordinates  9,  ^  the  flux  coordinate,  and 
X  the  magnetic  potential.   Consider  an  electrostatic  per- 
turbation j6  ^  j6{X  >$)  e   e    where  co  is  taken  to  be  small 
compared  to  the  frequencies  characteristic  of  the  guiding 
center  motion,  i.e.,  the  bounce  frequency  for  trapped 
particles  and  the  frequency  around  the  toroidal  minor  axis 
for  untrapped  particles.   In  this  limit  the  particle  sees 
only  the  time  average  of  the  potential  over  its  orbit,  and 
for  a  Maxwellian  equilibrium  the  perturbed  distribution  is 
given  by: 

,  mT  1  ^^o 


CD  +  mv -j 

where  v,  is  the  usual  guiding  center  drift  (d0/dt)  and  the 
bars  indicate  a  time  average  over  the  unperturbed  orbit. 
For  untrapped  particles  in  the  presence  of  shear  j6   vanishes, 

Using  quasineutrality  we  may  derive  a  variational 

2 
expression  which  is  an  extremal  for  tu  in  the  cases 

—  -J—  »  ^  -^  or  T.  =  T  .   A  necessary  and  sufficient 

condition  for  stability  is  then  that  v,  must  be  favorable 

for  all  particles,  or  all  trapped  particles  in  the  presence 
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of  shear.  I.e.,  ■^)    ^  <  o.   Growth  rates  are  of  order 

E  ,^- 
a^/r  smaller  than  hydrodynamlc .   Conventional  multipoles 

violate  this  condition  for  barely  trapped  particles. 

However,  we  show  that  a  proper  minimum  J  geometry  may  be 

constructed  by  modification  of  the  multipole  field. 
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Guiding  Center  Stability 
D.  B.  NELSON 
Oak  Ridge  National  Laboratory 

Using  a  variational  principle  for  the  guiding  center 
plasma  we  examined  the  conjectured  theorem  that  at  suitably 
low  p,  interchange  stability  plus  local  stability  is 
sufficient  for  absolute  stability.   The  geometry  chosen 
for  this  examination  is  a  mirror  machine  with  insulated 
ends,  so  that  the  magnetic  lines  are  free  rather  than  tied 
at  the  ends. 

Although  we  have  not  yet  proven  this  theorem  in  full 
generality  we  have  shown  that  for  the  subclass  of  equilibria 
which,  roughly  speaking,  admit  no  neutral  interchanges  the 
theorem  is  true  if  the  magnetic  field  and  plasma  pressure 
are  smooth  enough.   The  requirement  of  low  p  results  from 
integrating  the  small  pressure  derivative  in  from  the  plasma 
boundaries.   More  exactly,  if 

mm  — A^    ^   >  0 
U   /U^dv 

where  5  ^  is  the  second  variation  and  the  minimum  is  taken 
over  all  interchanges,  then  the  equilibrium  is  stable  if 

o 

it  is  locally  stable  and  if  V(pj_  +  B  /2 )  and  the  tensors  VB 

2        2 
and  5  (pj_  +  B  /2)/3x.Sx.  are  sufficiently  small. 
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The  method  of  proof  uses  the  fact  that  any  variation 
U  can  be  decomposed  Into  an  Interchange  U"^  plus  a  varia- 
tion U  which  vanishes  on  one  end.   Then  the  second 
variation,  as  a  quadratic  form,  becomes 

5^?(U)  =:  Q(U°)  +  Q(U°,U^)  +  Q(U^) 

Now  Q(U  )  is  positive  by  hypothesis.   We  have  shown  that 
Q(U  )  is  positive,  proving  a  theorem  that  for  a  smooth 
enough  (or  short  enough)  plasma  with  one  end  tied,  local 
stability  suffices  for  absolute  stability.   The  work 
remaining  to  complete  the  proof  involves  bounding  the 
cross  terms  Q(U°,U"^).   We  have  already  done  this  for  the 
case  given  above. 
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Absolute  Loss-Cone  Instabilities* 
R.  A.  DORY  and  G.  E.  GUEST 
Oak  Ridge  National  Laboratory 

Rosenbluth  and  Post   first  analyzed  the  high-frequency 
electrostatic  instabilities  resulting  from  Inverted  distri- 
bution of  gyration  energy  caused  by  the  loss-cone  in 
mirror  traps  or  by  velocity-dependent  loss  mechanisms  in 
other  configurations.   Because  the  modes  were  found  to  be 
convective  in  nature,  it  was  possible  to  find  heuristic 
sets  of  mirror  trap  parameters  that  provide  adequate 
stability  for  fusion  purposes.   The  inverted  distributions 
also  cause  other  instabilities  which  occur  at  lower  density 
and  which  are  associated  with  the  gyration  frequency  of  the 

Ions.   These  were  investigated  earlier  by  the  present 

2 
authors,   but  the  propagation  properties  of  the  waves  were 

not  determined.   We  have  extended  the  calculations  and 

have  found  that  under  a  range  of  circumstances,  these 

modes  are  absolute  rather  than  convective,  so  that  the 

stabilization  Invoked  for  the  Rosenbluth-Post  modes  does 

not  apply  here.   We  present  a  preliminary  investigation  of 

the  parameter  ranges  for  which  the  instabilities  are 

absolute.   It  is  found  that  even  very  smooth  inverted 

distributions  lead  to  absolute  instabilities.   To  obtain 

a  degree  of  stabilization  it  appears  likely  that  finite 
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length  effects  will  be  needed;  either  through  restric- 
tions on  machine  length  so  that  wavelengths  are  bounded 
and  Landau  damping  will  be  effective  (electron  heating 
may  be  necessary  to  provide  sufficient  damping),  or  If 
these  length  requirements  are  too  stringent,  It  may 
become  necessary  to  treat  the  finite  boundary  value 
problems  In  detail  In  order  to  find  new  stabilizing 
mechanisms . 


Research  sponsored  by  the  U.  S.  Atomic  Energy 
Commission  under  contract  with  the  Union  Carbide  Corporation. 

M.  N.  Rosenbluth  and  R.  P.  Post,  Phys.  Fluids  8, 
5^7  (1965). 

2 

G.  E.  Guest  and  R.  A.  Dory,  Phys.  Fluids  8,  I853 

(1965). 
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Flute-Like  "Loss-Cone"  Instabilities  In 
Multl -Component  Plasmas* 
W.  M.  FARR  and  R.  L.  BUDWINE 
Oak  Ridge  National  Laboratory 

Detailed  calculations  have  been  made  to  determine 
Instability  thresholds  for  high-frequency  flute-like 
oscillations  In  plasmas  containing  two  distinct  groups 
of  ions,  one  of  which  is  relatively  cold  and  well 
thermallzed  and  the  other  of  which  is  of  the  class  of 
"loss-cone"  distributions  studied  by  Guest  and  Dory. 

The  results  are  displayed  in  the  parameter  space  gener- 

2    2 
ated  by  N^Ntt  and  o)  /co   ,  where  N^  and   N^^  are  the  cold- 

and  hot-ion  densities,  cd  is  the  plasma  frequency 
corresponding  to  the  total  ion  density,  and  cd  is  the 
ion  gyrofrequency .   These  stability  boundaries,  corres- 
ponding to  the  "most  unstable"  ranges  of  wavelength,  are 

compared  to  the  less  specific  thresholds  given  earlier 

2 
by  Perlstein  et  al . 


Research  sponsored  by  the  U.  S.  Atomic  Energy 
Commission  under  contract  with  the  Union  Carbide  Corporation, 

"'"  G.  E.  Guest  and  R.  A.  Dory,  Phys.  Fluids  8^,  I853 
(1965). 

^  L.  D.  Perlstein,  M.  N.  Rosenbluth  and  D.  B.  Chang, 
Phys.  Fluids  9,  953  (1966). 
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Scattering  from  Magnetic  Fluctuations"**" 

WILLIAM  B.  THOMPSON 

University  of  California  at  San  Diego 
La  Jolla,  California 

Scattering  of  optical  and  radio  frequency  radiation 
from  density  fluctuations  In  a  plasma  has  been  extensively 
studied  and  bids  fair  to  become  a  useful  diagnostic  tool. 
Recent  observations  of  Faraday  rotation  in  the  optical 
region  imply  that  fluctuations  in  the  magnetic  field 
should  give  rise  to  a  further  scattering  which  might  also 
prove  useful  in  the  study  of  turbulent  plasmas. 

The  usual  scattering  results  may  be  obtained  from  the 
expression  for  the  scattered  spectrum  due  to  a  distribution 
of  polarization  P  over  a  volume  V  distance  R  from  the 
point  of  observation 

I(a))da3  =  ^-|   d^r'<^P^(r,CD)  .p^(o,0D)^e°  ^  ^ 

where  Pj_  is  the  projection  of  P  normal  to  R. 

One  writes  P  =  -(cd  /od  )E  and  oo  =  — —  [n  +5n(x-t)] 

^p'-'o  p  m'-o  ^^     '  ^ 

whence   follows 

-^ do^  =  n^Vg   87r[e  /mc    ]         , 

o  L  _ 

<(8n/n)2>  f(cD-a)^,k-k^).(Ejk^,a)^)    -  E^.R)^ 
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where   <6n^)f(cu,k)  =  j  <6n(r  .a3)6n-^(o,cD)  \  e"^-^"  ~  d  r'dcu 
so  that  f  Is  a  normalized  correlation  function. 

In  a  magnetic  field,  B,  the  polatlzatlon  may  be 
written 

P=  -0.^0.2  [E  +  gExb+  0(g)2] 

and  In  the  presence  of  fluctuations  as 

5P  =.  -  a.%^j  ^  [E  +  ^  E  X  b]  +  ^  [^  E  X  b  +  E  X  5b  ] 

p'     L  n    ^;     CD  -     -  -      CO    B   - 

eB 
where  fi  Is  the  gyro  frequency  In  the  field  B,  —  and  b, 

the  unit  vector  along  B.   Of  the  terms  In  here,  the  second 

>^ 

gives  a  depolarized  modification  ~  (fi/co)   the  usual 
scattered  wave,  and  Is  perhaps  negligible,  and  the  second 
and  third,  which  arise  from  magnetic  fluctuations,  are 
orthogonal  to  the  first.   This  orthogonality  persists  in 
the  scattered  field  only  in  special  directions  (e.g., 
backward,  forward.')  of  scattering.   In  these  special  direc- 
tions the  field  fluctuation  effect  may  be  separated  from 
the  rest  by  controlling  the  polarization  of  the  observed 
wave.   The  added  scattered  intensities  have  the  form  (for 
back  scattering) 

I^I^  =  nVg  87r(e%c2)^<6B^VB^  g (00-03^, k-k^) 
Ig/I^  =  nVg  87r(e%c2)2<502>  h(co-a3^,k-k^) 


where  cd,  k,  cd  ,  k  =  scattered  and  Incident  frequency 
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and  wave  number 

(g,h)    =    (fi/cD)^n(I,h) 

and 

,3 


r 


<6B^>g   =    f  <6B(r,a))5B(o,u))^   e^-*-   d' 

<60^>h   =    J  <(e^  X    5h(r,a))-EQ  x    6h(o,cu)      e^-"~    d\ 

g  and  h  have  order  of  magnitude  (n/tu)  ni"^  where  Ji   is  the 
correlation  length  and,  although  the  frequency  ratio  Is 
necessarily  small,  the  factor  n£     may  be  very  large.   Thus 
magnetic  scattering  may  be  as  large  as,  or  much  larger 
than  scattering  from  density  fluctuations.   It  can  be 
shown  to  be  small  in  thermal  equilibrium,  and  is  probably 
small  in  low  p  plasmas;  however,  in  high  p  turbulent 
plasmas  it  may  be  of  considerable  importance. 


Work  supported  in  part  by  the  U.  S.  Atomic  Energy 
Commission. 
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Bremsstrahlung  Emission  in  Plasmas 
ALAN  OPPENHEIM 
Polytechnic  Institute  of  Brooklyn 

We  develop  a  formula  for  the  power  density  of 
Bremsstrahlung  radiation  from  binary  interactions 
emitted  from  a  fully  ionized  non-relativistic  plasma. 
The  calculation  is  strictly  classical.   The  resulting 
formula  contains  a  linear  variation  with  temperature 
rather  than  the  variation  with  the  square  root  of  the 
temperature  as  computed  from  Kramer's  results.   The  two 
formulas  are  numerically  equal  at  4700  K.   Above  this 
temperature  the  linear  formula  clearly  predicts  a  higher 
rate  of  radiation  emission. 
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Boundary  Value  Problem  for  Bernstein  Resonances 

HAROLD  WEITZNER 

New  York  University 
Courant  Institute  of  Mathematical  Sciences 

A  plasma  varying  with  only  one  spatial  coordinate 
described  by  the  linearized  Vlasov  equation  with  an 
unperturbed  magnetic  field  perpendicular  to  the  direction 
of  spatial  variation  is  considered.   The  electric  field, 
with  only  the  longitudinal  part  Included,  is  given  by 
Polsson's  equation.   A  method  of  solution  is  described 
for  the  mixed  initial  value -boundary  value  problem  for  the 
case  of  a  plasma  bounded  by  one  wall  only.   When  the  plasma 
is  driven  by  a  harmonically  oscillating  electric  field  at 
the  wall,  the  behavior  of  the  electric  field  for  large 
time  is  examined.   The  initial  data  also  contribute  to 
the  long-time  limit  of  the  electric  field,  and  this  part 
is  given.   The  limiting  value  of  the  electric  field  for 
large  separation  from  the  wall  is  described  as  a  function 
of  driving  frequency.   The  line  shape  as  a  function  of 
frequency  is  given  near  resonance,  and  the  full  time 
behavior  at  resonance  is  examined. 
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Echo  Phenomena  Associated  with  Landau  Damping 

T.  M.  O'NEIL 

General  Atomic  Division 
General  Dynamics  Corporation 

R.  W.  GOULD 

California  Institute  of  Technology 


A  calculation  Is  presented  which  predicts  the  fol- 
lowing echo  phenomena.   A  wave  of  wave  length  k..  Is 
excited  In  a  plasma  and  then  Landau  damps  away;  after  a 
time  T,  a  wave  of  wave  length  kp  Is  excited  In  the  plasma 
and  It  also  Landau  damps  away;  then  after  a  time 
t  =  T[kp/ (kp-k-,  )  ]  ,  a  third  wave  (i.e.,  the  echo)  spon- 
taneously appears  In  the  plasma. 

The  basic  mechanism  behind  these  phenomena  can  be 
easily  understood.   When  the  first  wave  damps  away,  it 
leaves  a  perturbation  in  the  particle  distribution  function 
of  the  form  f,  (v)cos(k-,x  -  k-,vt).   There  is  no  electric 
field  associated  with  this  perturbation,  since  for  large 
t  the  perturbation  Is  a  rapidly  oscillating  function  of  v, 
and  a  velocity  Integral  over  the  perturbation  (i.e.,  the 
charge  density)  will  phase  mix  to  zero.   The  second  wave 
modulates  the  unperturbed  particle  distribution  and  leaves 
a  similar  perturbation,  but  it  also  modulates  the  perturba- 
tion due  to  the  first  wave,  leaving  a  second  order  perturba- 
tion of  the  form  f,  (v)f,  (v)cos[k,x  -  k-,  vt  ]cos[k^x  -k^v(t-T)] 
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This  perturbation  can  be  rewritten  as 
^k  ^^^^k  ^^^  cos[(k^+k2)x  -  (k-^  +  k2)vt  +  k^vl]   + 
cos[(k^-k2)x  +  (k^  -k^)vt  -  k^vT  ]  .   In  the  second 
cosine,  the  coefficient  of  v  will  vanish  when 
t  =  T[k2(k2  -k^)];  so  at  this  time  the  second  term  will 
no  longer  be  a  rapidly  oscillating  function  of  v.   Con- 
sequently, a  velocity  Integral  over  this  term  will  not 
phase  mix  to  zero,  and  an  electric  field  will  reappear 
In  the  plasma.   Of  course,  this  qualitative  explanation 
has  been  verified  by  a  perturbation  solution  of  the  Vlasov 
equation. 


This  research  was  sponsored  by  the  Defense  Support 
Agency  under  contract  DA  49-146-XZ-486. 
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Computation  of  Equilibrium  Configurations 

of  a  Toroidal  Plasma 

S.  FISHER  and  J.  KILLEEN 

Lawrence  Radiation  Laboratory,  University  of  California 

Livermore,  California 

We  have  developed  a  Fortran  program  for  computing  static 
hydromagnetic  equilibria  for  a  toroidal  plasma  with  scalar 
pressure.   We  assume  that  the  system  is  axially  symmetric  so 
we  can  apply  these  calculations  to  such  configurations  as 
the  levitron  or  floating  ring  devices.   Grad  and  Rubin 
have  shown  that  a  stream  function  i'{r,z)   can  be  introduced 
which  satisfies  the  equation 

OZ 


Where  f(T/')  =  p  and  g{'f)   =  ^B^'  i.e.,  the  pressure  and  azi- 

muthal  component  of  the  field  are  functions  of  f   only.   This 

2 
equation  was  solved  analytically  for  a  special  form  of  f 

and  g.   We  solve  the  above  equation  by  finite-difference 

methods  using  an  alternating  direction  implicit  method 

developed  for  an  earlier  problem  with  the  same  equation. 

The  choice  of  f  and  g  is  quite  flexible  and  various  profiles 

can  be  studied.   We  follow  the  procedure  of  the  earlier 

3 

problem,   and  let  i/  =  i/     +  if    ,   where  if     is  the  stream  func- 

c    p         c 

tion  of  the  vacuum  magnetic  field  and  i      is  the  stream  func- 

IT 

tion  of  the  magnetic  field  due  to  the  plasma.   The  function 
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if    (r,z)  Is  then  given  by  the  experimental  configuration 

and  Is  computed  from  the  coil  elements.    We  choose  fields 

with  closed  constant-''/'   surfaces  and  assume  that  the  plasma 

c 

is  contained  between  two  such  surfaces.   The  solution  for 
i/    (r,z)  is  by  iteration  of  the  difference  equation,  and 
a  succession  of  equilibrium  solutions  can  be  found  for 
increasing  values  of  the  pressure  maximum. 


Work  performed  under  the  auspices  of  the  U.S.  Atomic 
Energy  Commission. 

H.  Grad  and  H.  Rubin,  Proceedings  of  the  Second 
International  Conference  on  the  Peaceful  Uses  of  Atomic 
Energy,  31 »  190  (1958),  United  Nations,  N.Y. 

E.  W.  Lalng,  S.  J.  Roberts  and  R.  T.  P.  Whipple, 
J.  Nucl.  Energy,  Part  C  1,   49  (1959)- 

^  J.  Killeen  and  K.J.  Whlteman,  Phys .  Fluids  9,  l846 
(1966). 

u 

J.  Killeen,  H.  P.  Purth,  R.  P.  Prels,  "Calculation 

of  Toroidal  Magnetic  Field  Configurations",  UCRL-5OI6I, 

January  1967- 
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Toroidal  Magnetic  Fields  with  a  Circular  Magnetic  Axis 

GEORG  KNORR 

Department  of  Physics 
University  of  California,  Los  Angeles 

ECKHARD  REBHAN 

Institute  for  Plasma  Physics,  Garchlng,  Germany 

The  general  solution  of  Laplace's  equation  for  vacuum 
magnetic  fields  generated  by  arbitrary  surface  currents  on 
a  circular  torold  Is  given  In  toroidal  coordinates.   Mag- 
netic surfaces  are  constructed  around  a  circular  axis  as 
an  expansion  in  powers  of  the  aspect  ratio.   For  a  simple 
model  this  expansion  can  be  extended  in  principle  to 
arbitrarily  high  orders.   Whenever  the  rotational  transform 
t-/2-w   on  the  axis  assumes  rational  values  ra/n  (m  <  n  and 
without  common  divider)  there  might,  but  must  not,  appear 
resonances  (e.g.,  by  stray  fields  with  periodicity  n)  which 
lead  to  a  destruction  of  the  surfaces.   The  values  of  v" 
and  6/27r  are  computed  on  the  axis  and  expressed  by  parameters 
of  the  field.   On  the  other  hand,  the  field  generating  sur- 
face currents  can  also  be  expressed  by  the  field  parameters 
thus  relating  favorable  field  properties  directly  to  the 
currents.   It  is  shown  that  fields  with  arbitrarily  deep 
magnetic  wells  and  large  L   (and  shear)  can  be  produced  by 
relatively  simple  helical  coils. 
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Toroidal  Magnetic  Well  Configurations 

HAROLD  GRAD 

New  York  University 
Courant  Institute  of  Mathematical  Sciences 


We  analyze  the  following  statement: 
"A  plasma  Is  dlamagnetic.   It  Is  therefore  stable 
In  a  magnetic  well.   A  magnetic  well  Is  topologl- 
cally  Impossible  In  a  torus.   We  must  therefore 
settle  for  an  average  well  to  obtain  stability 
In  a  torus." 

1.  Minimum  B  (as  It  Is  Interpreted  In  a  mirror 
machine)  Is  not  related  to  plasma  stability. 

2.  Minimum  average  B  (/  di/B,  V",  etc.)  has  almost 
no  relation  to  stability. 

3.  Minimum  average  B  is  not  a  weakened  form  of 
minimum  B;  they  are  essentially  unrelated. 

4.  There  are  simple  examples  of  a  plasma  which  Is 
paramagnetic .   Whether  it  is  paramagnetic  or  dlamagnetic 
has  no  relation  to  its  stability  either  in  a  well  or  not. 

5-   There  are  very  simple  toroidal  vacuum  magnetic 
wells.   When  combined  with  a  proper  distribution  function 
[Taylor's  f{e,\i)]    they  are  stable  (5W  >  O).   No  considera- 
tions of  shear,  line  curvature,  or  V"  enter. 
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Equilibrium  Electric  Fields  In  Multlpoles 
T.  K.  FOWLER 
General  Atomic  Division/General  Dynamics  Corporation 

Low-p  equilibrium  configurations  in  symmetric  toroidal 
magnetic  geometries,  such  as  multlpoles,  are  examined  when 
toroidal  effects  are  Important.   In  this  regime,  locally 
Maxwellian  equilibrium  distributions  require  an  equilibrium 
ele'ctrlc  field.   For  T  «  T.,  the  extreme  cases  are: 

(1)  ions  held  electrostatically  with  potential  ep  =   .T^ ,  and 

(2)  ions  held  magnetically  with  potential  ep   =  ^^i/'^p^M^'^e 
where  R.  is  the  ion  gyroradius  and  R  and  R^^  are  the  minor 
and  major  radii.   The  first  case  is  characterized  by  zero  ^ 
ion  current,  small  total  angular  momentum,  uniform  pressure 
along  magnetic  lines  and  E-B  =  0.   In  the  second  case,  ions 
carry  current  and  a  sizeable  momentum,  and  there  must  be  a 
small  electric  field  along  B  which  pushes  electrons  toward 
outer  radii  where  ions  are  driven  by  centrifugal  forces  (the 
toroidal  effect).   Some  implications  of  these  results  will 
be  discussed. 


-^ 

This  work  was  carried  out  under  a  joint  General  Atomic- 
Texas  Atomic  Energy  Research  Foundation  program  on  controlled 
thermonuclear  research. 
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Closed  Line  Toroidal  Magnetic  Field  Configurations 

H.  P.  FURTH"^  and  J.  KILLEEN 

Lawrence  Radiation  Laboratory,  University  of  California 

Livermore ,  California 


There  is  some  Interest  in  toroidal  magnetic  field  con- 
figurations in  which  the  field  lines  close  upon  themselves 
and  have  a  maximum  of  ^  di/B  in  the  plasma  region.   Floating 
ring  devices  have  this  property.   A  linear  periodic  system 
with  this  property,  but  without  a  floating  conductor  was 

discussed  by  Furth  and  Rosenbluth.    Recently  a  toroidal 

2 
configuration  with  this  property  was  studied  by  Taylor.    We 

have  also  been  studying  toroidal  configurations  without 

floating  conductors.   To  do  this  we  have  found  it  useful  to 

use  a  spherical  coordinate  system  p,  6,   j6.      The  magnetic 

field  consists  of  two  parts,  the  ordinary  toroidal  field 

given  by  B^  =  ? — n   and  a  field  B  =  VX  ,  where  the  function 

°      "^   jS   p  sxn  9 

y(p,9,j25)  is  a  sum  of  various  spherical  harmonics.   Approximate 
analytical  work  has  indicated  that  a  stable  §   d^B  region  can 
be  obtained.   To  study  this  exactly  we  have  developed  a  new 
Fortran  program  called  TUBE  V  which  solves  for  the  field  line 
trajectories  in  a  torus,  but  in  spherical  coordinates.   This 
technique  of  employing  spherical  harmonics  is  useful  in 
studying  more  general  stellerator  fields.   We  shall  present 
results  of  some  closed  line  configurations. 
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Work  performed  under  the  auspices  of  the  U.S.  Atomic 
Energy  Commission. 

"^Now  at  Plasma  Physics  Laboratory,  Princeton  University, 
p.  0.  Box  451,  Princeton,  New  Jersey. 

■^H.  P.  Furth  and  M.  N.  Rosenbluth,  Phys.  Fluids  7,  764 
(1964). 

J.  B.  Taylor,  "Anti-Synmietric  Toroidal  Containment 
Systems,"  Preprint  Culham  Laboratory,  1 
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Diffusion  in  a  Cylindrical  Force  Free  Field 

M.  BINEAU 
New  York  University 
Courant  Institute  of  Mathematical  Sciences 

The  choice  of  an  arbitrary  function  f(r)  in  an 

Internal  (0,R)  or  (R-j^,!^)  as  the  ratio  of  the  current 

density  to  the  magnetic  field  determines  a  unique 

cylindrically  symmetric  force  free  field  in  various 

situations.   This  occurs  for  a  small  e,  where  e  measures 

the  ratio  of  the  magnetic  field  induced  by  the  current 

to  the  given  vacuum  field.   This  occurs  also  for 

arbitrary  e  and  appropriate  boundary  conditions  leading 

to  a  Vol terra  equation.   The  medium  supporting  the  force 

free  field  is  assumed  to  be  resistive  and  its  behavior 

governed  by  Ohm's  law,  the  force  free  condition  and 

Maxwell's  equations.   When  the  externally  Induced  part 

of  the  field  is  kept  constant,  the  velocity  v  of  the 

medium,  which  is  also  the  drift  velocity  of  individual 

particles  v  =  — ^^  ,  is  obtained  explicitly  after 

B 

solving  equations  in  integral  form  for  the  time 
derivative  of  the  current  density.   For  small  s  the 
solution  can  be  obtained  as  an  expansion  In  powers  of  e. 
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The  lowest  order  term  in  e  is 
R  X 

V  =  E^r  J     ^^   (f(x)  -  \    J  t   f(t)  dt)  dx  for 

r  ^0 

a  force  free  field  In  a  domain  (0,R)  and  resistivity 

T].   For  a  force  free  field  In  a  domain  (R-,,Rp)  with 

axis  Oz  (hard  core  conf Iguraglon) ,  two  natural  boundary 

aB^(R^) 
conditions  can  be  considered,  namely  -vr —  =  0  or 

SB  (  R2 ) 

— -rr =  0.   The  corresponding  velocities  are  to 

lowest  order  In  e, 

V  =  ,  ^  /  ££(il^  and  V  =  -n  ^  /■"'  S£i^^  . 

Bz  Bz 


These  formulae  may  be  compared  with  the  relative  velocity 
V  -  v^  =  -rif  tg—  ,  where  v..  Is  the  velocity  -rr-  which 

M         xj  M  d  t 

Z 

preserves  constant  flux  In  a  circle  of  radius  R.   Note 
that  the  velocity  is  going  into  the  domain  if  the 
current  in  the  medium  flows  in  the  opposite  direction 
to  the  current  in  the  hard  core. 
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Energy  Variational  Principles  for  Steady 
Hydromagnetic  Plows  in  a  Toroid 
M.  D.  KRUSKAL  and  C  H.  SU 
Plasma  Physics  Laboratory,  Princeton  University 

We  consider  ideal  hydromagnetic  flows  in  a  toroid. 
The  work  is  an  extension  of  Kruskal  and  Kulsrud's  [l] 
study  on  static  equilibria.   It  is  shown  here  that  a 
hydromagnetic  flow  is  steady  if  and  only  if  it  makes  the 
total  energy  of  the  system  variationally  stationary  under 
the  constraint  that  the  five  surface  functions  not  vary: 

s(m)  ,  i^(m)  ,    (m)  ,  W(m)  ,  N(m) 

representing  specific  entropy,  long-  and  short -way 
magnetic  fluxes,  and  two  constants  of  motion  pertinent 
to  flow  velocity.   Here  m  =  M(x)  is  the  magnetic  surface 
enclosing  mass  ra.   A  detailed  discussion  on  the  last  two 
constants  of  motion  W(m)  and  N(m)  is  given  based  on  recent 
work  of  Newcomb  [2].   Two  distinct  energy  variational 
principles  are  constructed.   Finally  it  is  shown  that 
one  can  derive  from  our  energy  variational  principle  the 
Lagranglan  variational  principle  of  Greene  and  Karlson  [3], 


M.D.  Kruskal  and  R.M.  Kulsrud,  Phys .  Fluids  1,  265 
(I950).  — 

2 
W-A.  Newcomb,  UCRL-14884,  University  of  California, 
Llvermore,  California  (1966) . 
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^J.  M.  Greene  and  E.  T.  Karlson,  Princeton  Plasma 
Physics  Laboratory  Report  MATT-478  (I966). 
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Acceleration  of  Macropartlcles  to  Very  High 

Velocities  by  Megagauss  Fields 

F.  WINTERBERG 

Desert  Research  Institute  and  Department  of  Physics 
University  of  Nevada,  Reno,  Nevada 

A  new,  very  promising  method  for  the  acceleration 
of  macropartlcles  to  extremely  high  velocities  by  mega- 
gauss fields  Is  proposed.   By  this  method.  It  seems  to 
be  possible  to  reach  velocities  of  the  order  of  lO'  cm/sec 

o 

and  perhaps  10  cm/sec.   It  was  previously  shown  by  the 
author  that  a  small  thermonuclear  explosion  can  be 
triggered  by  the  Impact  of  macropartlcles  with  velocities 
exceeding  10'  cm/sec  in  dense  T-D.   The  described  new 
method  of  macroparticle  acceleration  may  therefore  lead 
to  the  controlled  release  of  thermonuclear  energy.   This 
method  of  acceleration  is  based  on  a  new  and  very  promising 
method  for  the  generation  of  magagauss  fields.   In  contrast 
to  previous  methods  using  high  explosives,  the  new  method 
to  generate  megagauss  fields  employs  a  fast  moving  pro- 
jectile for  flux  compression.   The  projectile  is  fired 
either  by  an  ordinary  or  ultimately  by  a  light  gas  gun. 


Case  Institute  of  Technology,  Plasma  Research 
Program,  Technical  Report  A21,  June  I963  and  Z.  Naturf. 
19a,  231  (1964). 
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On-Line  Solution  of  Loss-Cone  Dispersion  Relations 

BURTON  D.  FRIED 
University  of  California,  Los  Angeles  and  TRW  Systems 

JAMES  E.  McCUNE 
Massachusetts  Institute  of  Technology 

A  convenient  technique  for  solving  dispersion  equations 
of  the  kind  encountered  in  plasma  physics  has  been  described 
previously.    Application  of  this  method  to  the  complicated 
dispersion  equations  associated  with  a  collisionless,  mag- 
netized plasma  is  greatly  facilitated  by  an  on-line  digital 
computer  which  provides  easy  access  to  functional  operations 

In  the  complex  plane,  as  well  as  immediate  displays  of  arcs 

2 
or  contours.    New  features  and  recent  refinements  of  this 

method  are  discussed  and  the  technique  is  then  applied  to  the 

case  of  the  Rosenbluth-Post  "loss  cone"  instability  in  its 

most  virulent  form,  i.e.,  when  the  distribution  function  is 

2   2 
of  the  form  exp ( -v  /a^)6(vj^  -  aj^).   The  transition  from  the 

case  of  small  x  =  (kj_a^(jo  .  ) ,  where  a  few  terms  of  the  Bessel 

function  series  dominate  the  behavior,  to  the  case  of  large  x 

is  illustrated.   The  relation  of  these  exact  solutions  to 

those  obtained  in  the  high  density  limiting  case  studied  by 

Rosenbluth  and  Post   is  discussed. 


J.E.  McCune,  Phys.  Fluids  9,    2082  (1966)  and  B.D.  Pried, 
"On-Line  Root  Finding  in  the  Complex  Plane,"  Report  No. 
9990-7308-ROOOO  TRW  Systems,  Redondo  Beach,  Calif. 
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B.D.  Fried,  "Solving 
of  "On-line  Computing",  ed. 
New  York,  I967.   Also  G.J. 
On-line  System,"  Appendix, 


Mathematical  Problems,"  Chap.  VI 

by  W.  Karplus,  McGraw  Hill, 
Culler,  "Users  Manual  for  an 
loc .  cit . 


3 

"^M.N.  Rosenbluth  and  R. 

(1965)  and  9,    730  (1966). 


F.  Post,  Phys  Fluids  8,  547 
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Computer  Simulation  of  Toroidal  Plasma  Confinement 
Including  Dlsslpative  and  Inertlal  Effects 
N.  K.  WINSOR,  J.  L.  JOHNSON,"^  and  J.  M.  DAWSON 
Plasma  Physics  Laboratory,  Princeton  University 

A  computer  program  Is  being  developed  to  study  the 
confinement  of  low-P  plasmas  from  the  fluid  point  of 
view  In  simple  axlally  symmetric  configurations.   It 
should  be  possible  to  study  ideal  stellarators  (the  model 
of  Knorr  [1]  and  Karlson  [2]),  toroidal  pinches,  and  the 
levitron  on  the  present  model.   With  minor  modifications 
it  will  be  possible  to  study  spherators  and  multipoles. 
The  fluid  equations  adapted  may  include  resistivity, 
viscous  and  gyro-viscous  effects  (general  stress  tensor), 
the  Hall  term,  gradients  in  the  pressure  along  the  lines 
of  force,  and  the  effects  of  plasma  Inertia.   This  fluid 
approach  should  supplement  the  guiding  center  calculations 
of  Bishop  and  Smith  [3]  who  must  solve  a  self -consistent 
problem.   It  should  lead  to  a  better  understanding  of  the 
limitations  of  the  present  theories  of  plasma  equilibria 
and  the  plasma  losses  associated  with  classical  diffusion. 
It  should  help  to  clarify  the  role  played  by  the  mass 
motions  of  the  plasma  (streaming  along  field  lines  and 
plasma  rotation)  and  the  associated  dlsslpative  effects. 
The  method  for  attacking  this  problem  and  the  present  state 
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of  the  analysis  will  be  presented 


On  loan  from  Westinghouse  Research  Laboratories. 

"^  G.  Knorr,  Phys .  Fluids  3,  133^  (I965). 

2 

E.  T.  Karlson,  Princeton  Plasma  Physics  Laboratory 

Report  MATT -470  (I9b6). 

■^  A.  S.  Bishop  and  C  G.  Smith,  Phys.  Fluids  9,    1380 
(1966). 
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Nonlinear  Study  of  Vlasov's  Equation  for 
a  Special  Class  of  Distribution  Functions 
HERBERT  L.  BERK  and  KEITH  V.  ROBERTS 
University  of  California,  San  Diego 

Calculations  on  the  "water  bag  model"  of  a  two-stream 
Instability  have  been  made  by  a  new  numerical  method.   The 
distribution  function  Is  locally  randomized  but  large 
scale  nonlinear  waves  are  observed.   These  waves  are  due 
to  the  creation  of  holes  In  phase  space  that  formally 
behave  as  particles  of  negative  mass  since  holes  of  like 
charge  are  attractive.   On  a  longer  time  scale  the  system 
is  expected  to  approach  a  Permi  distribution. 
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Fourier  Hermlte  Expansion  of  the  Vlasov  Equation 

M.  R.  FEIX  and  F.  C.  GRANT 
The  College  of  William  and  Mary  and  NASA  Langley 

The  analytical  and  numerical  properties  of  the  double 
Fourier  Hermlte  expansion  of  the  Vlasov  equation  are  examined 
In  the  linear  and  the  nonlinear  case.   The  difficulties 
associated  with  the  closure  conditions  and  the  appearance 
of  the  mlcrostructure  In  velocity  space  are  shown  to  be 
rejected  at  larger  and  larger  time  where  N  (number  of  Hermlte 
polynomials  retained)  Increases.   Moreover  It  Is  shown  how  In 
the  linear  problem  the  two  limits  N  — ><»,  v  — >  0  with  Nv  >  1 
or  Nv  <  1  correspond  respectively  to  the  Landau  and  Van 
Kampen  results  when  we  introduce  a  small  but  finite  Fokker 
Planck  term.   The  Landau  limit  (Nv  >  1   N  — >  0°  v  — >  0)  Is 
much  more  Interesting  from  a  numerical  point  of  view  when 
nonlinear  numerical  calculations  are  considered. 

In  the  nonlinear  case  the  Fourier  Hermlte  expansion  is 
shown  to  correspond  to  a  double  quasi-linear,  long -wavelength 
expansion.   Numerical  results  are  presented  in  two  cases. 

(a)   Long  wavelengths  (kD  =  .05  and  harmonics),  a  =  .02 
(depth  of  the  initial  density  modulation).   The  small  thermal 
effect  destroys  the  exact  periodicity  of  the  cold  plasma  case, 
and  the  harmonics  (but  not  the  fundamental)  display  strong 
nonlinear  behavior.   Although  the  kD  are  very  small,  rather 
large  number  of  Hermlte  polynomials  (50)  are  necessary  to 


follow  the  system  during  long  time  (cd  t  =  80)  ,  indicating 
a  level  of  klneticity  much  stronger  in  the  nonlinear  case 
than  in  the  linear  one . 

(B)   Kinetic  regime  (kD  =    -5  and  harmonics),  the  Landau 
damping  is  found  to  decrease  when  the  amplitude  of  the 

field  is  increased  in  agreement  with  results  obtained  by 

-|  2 

Knorr  and  Armstrong  .   Numerical  results  indicate  that 

this  effect  is  a  quasi-linear  one  (i.e.,  does  not  involve 

the  first  harmonic).   A  discussion  of  the  role  of  this 

harmonic  is  presented  for  the  damping  case  showing  that 

although  wave-wave  interactions  are,  a  priori,  an  effect  of 

the  same  order  as  the  quasi-linear  one,  it  seems  reasonable 

to  disregard  it  for  a  small  damping  in  agreement  with  the 

numerical  results. 


■""  G.  Knorr,  Zeits.  fiir  Naturf orschung  l8a,  1304  (1963) 

2 

Armstrong,  APS  me< 

Boston,  1966,  paper  7P^ • 


Armstrong,  APS  meeting.  Division  of  Plasma  Physics, 
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Self-Consistent  Two -Dimensional  Equilibrium 

Calculations  Including  Particle  Loss 

J.  PREIDBERG  and  W.  GROSSMANN 

New  York  University 
Courant  Institute  of  Mathematical  Sciences 

We  discuss  computations  of  two-dimensional  equilibria 
applicable  for  both  low  and  high  p  mirror  devices.   Our 
aim  Is  to  study  the  effect  of  self  mirroring  due  to 
particle  currents.   Specifically,  we  are  Interested  In 
determining  the  self -consistent  particle  loss,  taking  into 
account  self  mirroring,  and  its  effect  on  equilibrium  flux 
surfaces.   Our  approach  consists  of  calculating  a  current 
from  an  equilibrium  solution  of  the  Vlasov  equation.   This 
current  is  then  used  with  Ampere's  law  to  compute  the  flux 
surfaces. 
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More  on  Toroidal  Geometry  with  Relatlvlstlc  Electron  Coils 

N.  C.  GHRISTOFILOS 

Lawrence  Radiation  Laboratory,  University  of  California 

Llvermore,  California 

The  combination  of  toroidal  geometry  with  true  magnetic 
well  created  by  relatlvlstlc  electron  colls  has  been  further 
studied.   The  required  strength  of  the  electron  colls  to 
satisfy  T.  K.  Fowler's  sufficient  condition  for  hydromagnetlc 
stability  has  been  previously  calculated. 

It  has  recently  been  observed  that  Fowler's  depth  of 

2   2 
the  well  (AB/b)  is  related  to  the  quantity  (cd^cd^  )  of  the 

relatlvlstlc  electron  coils  with  the  simple  Inequality 

(AB/B)  >  (0^0)2) 

Hence  the  sufficient  condition  for  hydromagnetlc  stability 

becomes 

T 
,  2/  2v    1    -^e 
(co^co^)  =  2  T  +  T. 


where  T  ,  T.  is  the  electron  and  ion  temperature  respectively 

of  the  plasma  confined  In  the  magnetic  well.   This  result 

2   2 
Indicates  that  the  value  of  (ow'to  )  is  independent  of  the 

shape  of  the  magnetic  well  created  by  the  relatlvlstlc 

2   2 
electron  coils.   For  T   ==  T.,  (^/<^  )  =  •25?  a  value  which 

has  been  achieved  in  the  Astron  experimentally.   At  this 

level  the  E-layer  is  still  free  of  instabilities.   Thus  a 

toroidal  geometry  of  closed  magnetic  lines  with  true  magnetic 
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well  created  by  a  large  liiimber  of  E-layers  disposed  in  a 
toroidal  geometry  is  now  in  the  realm  of  reality.   A  pre- 
liminary set  of  parameters  for  such  a  geometry  will  be 
presented . 


Work  performed  under  the  auspices  of  the  U.S.  Atomic 
Energy  Commission. 

"'"  H.  Christofilos,  The  Physics  of  Fluids  9,    1^25  (I966) 
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Interpretation  of  Experiments  on  Colllslonal  Drift  Modes 

B.  COPPI*  and  P.  PERKINS 
Plasma  Physics  Laboratory,  Princeton  University 


We  report  the  Interpretation  of  experiments  on 
low-temperature  alkali  plasmas  in  strong  magnetic  fields 
in  terms  of  colllslonal  drift  modes,  in  which  diffusion 
over  the  transverse  wavelength  resulting  from  ion-ion 
collisions,  plays  an  Important  role.   The  experimental 
results  agree  with  the  theoretical  predictions. 


On  leave  from  the  University  of  California  at 
San  Diego. 
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The  Inertlal  Confinement  of  Fusion  Gases: 

Theoretical  Aspects 

ROBERT  L.  HIRSCH 

ITT  Parnsworth  Research  Corporation, 
Port  Wayne,  Indiana 

A  theory  Is  presented  for  the  non-magnetic  Inertlal 
confinement  of  fusion  gases  In  spherical  geometry.   Assum- 
ing monoenergetlc  Ion  and  electron  distribution  functions, 
the  Polsson  equation  for  bipolar  charges  Is  solved  numeri- 
cally.  The  results  Indicate  that  within  a  broad  range  of 
currents  and  energies  periodic  solutions  are  obtainable. 
These  solutions  are  physically  Interpreted  as  the  alternate 
formation  of  virtual  anodes  and  virtual  cathodes  near  the 
center  of  a  hollow  spherical  cathode.   This  arrangement 
Involves  radial  ion  focusing,  thereby  creating  dense  regions 
of  high  kinetic  energy  ions  near  the  center  of  the  sphere. 
I.e.,  conditions  favorable  to  a  high  fusion  rate.   "Real 
gas"  modifications  to  the  theory  are  discussed.   In  particu- 
lar, a  mechanism  by  which  high  fusion  rates  enhance  ion 
trapping  is  presented.   A  means  of  direct  energy  conversion 
inherent  to  the  concept  is  also  described.   A  companion 
article  describes  encouraging  results  of  recent  experiments 
utilizing  this  concept. 
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Interaction  Between  Colllslonless  Shock  and 

Ion  Acoustic  Waves  in  a  Magnetic  Field 

N.  A.  KRALL  and  D.  L.  BOOK 

General  Atomic  Division  of  General  Dynamics  Corporation 
John  Jay  Hopkins  Laboratory  for  Pure  and  Applied  Science 

San  Diego,  California 

A  recently  proposed  model  to  explain  the  behavior  of 
turbulent  heating  experiments  by  means  of  a  spectrum  of 
ion  sound  waves  generated  by  diamagnetic  drifts  at  the 
shock  front  is  examined.   The  Vlasov  equation  is  used  to 
obtain  the  conditions  for  the  existence  of  instability, 
and  a  quasi -linear  equation  describing  the  effect  of  the 
unstable  spectrum  on  the  particle  distribution  functions 
is  derived.   Non-linear  corrections  are  studied  and  shown 
not  to  modify  the  picture  significantly.   It  is  concluded 
that  the  unstable  spectrum  is  not  built  up  to  a  level  high 
enough  to  provide  a  mechanism  for  the  particle  heating  and 
shock  form  observed. 


This  work  was  sponsored  by  the  Defense  Atomic 
Support  Agency  under  contract  DA-49-l46-XZ-53^ • 

"*"  R.  K.  Kurtmullaev,  Yu .  E.  Nesterikhin,  V.  I.  Pllsky 
and  R.  Z.  Sagdeev,  "Plasma  Heating  by  Colllslonless  Shock 
Waves,"  Plasma  Physics  and  Controlled  Nuclear  Fusion 
Research,  Vol.  II,  IAEA,  Vienna,  367-387  (1966). 


40 


The  Final  Solution  of  the  Problem  of 
Non-linear  Interactions 
R.  BUDWINe'''',  0.  C.  ELDRIDGE"^,  E.  G.  KARRIS'^", 
R.  SUGIHARA"*",  G.  M.  WALTERS"'""'" 

In  a  plasma  there  Is  a  finite,  and  usually  small, 
number  of  kinds  of  particles  ond  quasi -particles .   The 
Interaction  among  these  may  be  described  by  an  Interaction 
Hamlltonlan.   Once  this  is  known  the  problem  of  non-linear 
interactions  is  essentially  solved;  all  that  remains  is 
the  calculation  of  the  development  of  processes  of  physical 
interest . 

For  Instance,  if  Bo  =  0,  the  particles  are  electrons 
and  Ions,  and  the  quasi -particles  are  photons,  plasmons 
and  phonons.   The  coulomb  interaction  between  particles  and 
the  photon-particle  interactions  must  be  regarded  as  funda- 
mental and  the  others  (such  as  particle -plasmon,  plasmon- 
phonon,  etc.)  should  be  derived  from  them.   We  shall  discuss 
our  present  knowledge  of  these  Interactions  and  the  applica- 
tion to  derivation  of  the  quasi -linear  equations  and  higher 
approximations,  photon-photon  scattering,  radiation  from 
plasma  oscillations  and  acoustic  waves,  etc.,  etc. 

If  Bo  f   0,  there  are  several  additional  quasl-particles 
such  as  the  quanta  of  Alfven  waves,  whistlers,  etc.   We 
have  derived  the  partlcle-plasmon  and  plasmon-plasmon 
interaction  Hamiltonlans  and  used  them  in  the  derivation 
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of  quasi-linear  equations  and  the  wave-wave  coupling  -. 
corrections.   Also,  the  interaction  Hamiltonlan  between 
relatlvlstic  particles  and  (approximately)  transverse 
waves  has  been  derived  and  used  to  discuss  synchrotron 
radiation  and  the  cyclotron  maser. 


This  work  was  done  partly  at  the  Oak  Ridge  National 
Laboratory  (operated  by  the  Union  Carbide  Corporation  for 
the  U .  S.  Atomic  Energy  Commission)  and  the  University  of 
Tennessee  (under  contract  AT-(40-l )-2598  with  the  U.  S. 
Atomic  Energy  Commission). 

Oak  Ridge  National  Laboratory. 

University  of  Tennessee  and  Oak  Ridge  National 
Lab  orat  ory . 

''  University  of  Tennessee. 
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Non-Expandable ,  Non-Linear  Effects 

In  Weakly  Turbulent  Plasmas 

R.  E.  AAMODT 

Department  of  Physics,  University  of  Texas 

Austin,  Texas 

The  response  of  a  homogeneous  plasma  containing  a 
given  spectrum  of  longitudinal  waves  to  an  infinitesimal 
test  wave  is  found  in  a  local  approximation.   The  smear- 
ing of  the  wave-particle  resonance  and  the  resulting 
change  of  the  dispersion  relation  for  a  gentle  bump 
Instability  is  presented.   The  time  development  of  a 
weak  homogeneous  beam  in  this  system  is  examined  and  some 
new  non-linear  effects  important  in  beam  stabilization 
are  found . 
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Nonlinear  Interaction  of  Positive  and  Negative 
Energy  Modes  in  a  Plasma* 
B.  COPPl"*" 
Plasma  Physics  Laboratory,  Princeton  University 
M.  N.  ROSENBLUTH 
University  of  California  at  San  Diego 
R.  N.  SUDAN 
Plasma  Physics  Laboratory,  Princeton  University 


Nonlinear  instabilities  result  from  the  interaction 
of  positive  and  negative  energy  modes.   In  particular, 
a  high-temperature  plasma  confined  by  magnetic  mirrors  is 
considered.   For  frequencies  small  compared  to  the  ion 
gyrof requency,  in  the  linear  regime  two  flute  modes  are 
associated  with  the  magnetic  field  curvature  and  the 
density  gradient.   When  the  effect  of  the  ion-diamagnetic 
velocity  is  considered  both  these  modes  are  stabilized, 
but  one  is  a  positive  energy  mode  and  the  other  a  negative 
energy  mode  if  the  magnetic  curvature  is  unfavorable  to 
stability.   In  the  nonlinear  regime,  mode-mode  coupling 
Introduces  an  instability  which  upsets  this  finite 
Larmor  radius  stabilization.   We  also  consider  flute 


-  44  - 


modes  at  the  harmonics  of  the  ion  gyrof requency  deter- 
mined by  the  shape  of  the  distribution  function  and  the 
radial  density  gradient.   For  loss-cone  type  distribu- 
tions, we  observe  that  the  transition  from  positive  to 
negative  energy  modes  occurs  for  wavelengths  of  the  order 
of  the  ion  gyroradius.   For  Maxwellian  distributions  the 
same  transition  occurs  at  much  shorter  wavelengths. 
Therefore,  for  a  loss-cone  distribution  the  relevant 
nonlinear  instability  is  likely  to  cause  a  larger  particle 
diffusion  because  longer  wavelengths  are  involved.   From  a 
different  point  of  view,  this  instability  could  well 
explain  the   rf  bursts  at  harmonics  of  the  gyrof requency 
observed  in  several  mirror  experiments,  because  the  non- 
convectlve  nature  of  the  relevant  modes  may  allow  them  a 
sufficient  time  for  interaction. 


This  work  was  performed  in  part  under  the  auspices  of  the 
Air  Force  Office  of  Scientific  Research,  Contract  No. 
AP49(638)-155. 

'''On  leave  from  the  University  of  California  at  San  Diego. 


'On  1 


eave  from  Cornell  University,  Ithaca,  New  York 
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Development  of  a  Large  Amplitude  Wave  from  the 

Two-Stream  Instability 

JAMES  E.  DRUMMOND 

Boeing  Scientific  Research  Laboratories 
Seattle,  Washington 

The  large  amplitude  steady  state  waves  found  by 
Sen  are  re-calculated  and  found  to  result  from  distri- 
buted charges,  not  concentrated  dlpoles.   These  solutions 
are  then  perturbed  by  a  small  growth  rate  and  found  to 
develop  from  lower  amplitude  waves.   The  steady  state 
solutions  are  stable . 
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Quasi-Linear  Theory  for  Alfven  Waves 

FANG  C.  HOH 

Boeing  Scientific  Research  Laboratories 
Seattle,  Washington 


We  consider  the  following  equilibrium:   An  infinite 
homogeneous  plasma  is  immersed  in  an  uniforro  magnetic 
field   B  directed  along  the  z  axis.   The  plasma  is 
collision  free  and  has  arbitrary  pressure.   The  equilib- 
rium is  perturbed  by  a  plane  Alfven  wave  propagating 
along  the  z  axis  with  wavenumber  k.   The  Vlasov-Maxwell 
equations  are  used  and  are  Fourier  decomposed.   According 
to  the  quasi-linear  approximation,  we  only  keep  the  non- 
linear terms  in  the  diffusion  equation  that  describes  the 
time  change  of  the  Fourier  component  of  the  background 

distribution  function  f,     or  f  . 

k=o       o 

The  Independent  variables   v   and   v,  ,  the  particle 
velocities  in  the  direction  of  and  perpendicular  to  the 


magnetic  field  respectively  are  transformed  into  a  new 

^  ^-f  2v  V.  , 
z     z  A' 


2     P 
pair  of  variables;   ti  ^  v ,   and   ?=v,  +v-f2vv„. 


where   v^   is  the  Alfven  speed.   The  diffusion  equation 
then  reads : 
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where   D   is  a  diffusion  coefficient  containing  the 
square  of  the  wave  amplitude.   The  characteristic  line 
of  this  parabolic  differential  equation  is   ^  =  const, 
which  describes  an  elllps  in  the  {v^,    v    )    plane.   The 
unstable  particles  diffuse  along  the  contours  of  the 
ellipses  until   f    Is  constant  along  It,  i.e., 
hpyhf]   =   0. 

If  we  shoot  a  weak  beam  of  plasma  having  velocities 
close  to  V,  and  v  ^  (O  -  cD)/k  (O  =  ion  cyclotron  frequency 
and  CO  =  kv.  )  into  a  Maxwellian  plasma  (a  stable  equilib- 
rium), the  beam  ions  will  move  along  the  ^  =  const, 
contours  and  thereby  lose  their  Vj_  energy  which  partly 
goes  to  increase  their  v   energy  and  partly  to  excite  an 
Alfven  wave. 
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Expansion  of  a  Resistive  Plasmoid  in  a  Magnetic  Field* 

D.  K.  BHADRA 

General  Atomic  Division  of  General  Dynamics  Corporation 
John  Jay  Hopkins  Laboratory  for  Pure  and  Applied  Science 

San  Diego,  California 


Calculations  are  presented  which  show  the  effect  of 
resistivity  in  the  expansion  of  a  laser  produced  hot 
plasma  against  a  magnetic  field.   Without  resistivity, 
the  expanding  plasma  has  periodic  behavior  due  to 
bouncing  back  from  the  field;  but  with  resistivity,  there 
is  superimposed  a  slow  diffusion  (across  the  field)  on 
the  periodic  oscillations.   This  Is  illustrated  by 
considering  a  spherically  symmetric  plasma  produced  by 
a  giant  laser  pulse  (delivering  power  of  the  order  of 
10"^^W)  Illuminating  a  pertlcle  of  dimension  10   cm  and 
subsequently  following  the  time  development  of  its  radius, 
skin-depth  and  temperature  by  integrating  the  non-linear 
equations  of  motion  and  the  Ohm's  Law. 


This  research  was  sponsored  bv  the  Defense  Atomic  Support 
Agency  under  Contract  DA-49-146-XZ-4eb. 
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Kelvin-Helmholtz  Instability  In  Magneto-Plasmas 

GEORGE  SCHMIDT 
Stevens  Institute  of  Technology 

A  general  equation  describing  instabilities  in  a  low 
P  magneto  plasma  under  the  influence  of  nonuniform  electric 
field  and  gravitational  field,  both  perpendicular  to  B  , 
has  recently  been  derived  by  Stringer  and  Schmidt.    This 
eigenvalue  equation  that  contains  finite  gyro-radius 

effects  has  been  shown  to  reduce  in  the  appropriate  limits 

2 
into  the  equations  obtained  by  Rosenbluth  and  Simon  for 

high  density  and  low  density  plasmas.   It  has  also  been 
shown  that  in  the  very  high  density  and  zero  gyro-radius 
limit  the  equation  reduces  to  a  known  equation  in  hydro- 
dynamics that  describes  the  Kelvin-Helmholtz  instability. 
We  have  reduced   the  high  density  equation  by  using  a 
small  (but  not  zero)  gyro-radius  ordering  to  a  similar 
fluid  equation  and  obtained  necessary  conditions  for 
instability  in  terms  of  a  modified  Richardson  number.   In 
the  case  of  no  gravity  a  sufficient  condition  for  instability 
has  also  been  found.   Note  that  the  shear  flow  produced  by 
the  nonuniform  electric  field  may  induce  Instability  even 
in  the  presence  of  stabilizing  gravity. 


1.  T.  Stringer,  G.  Schmidt,  Plasma  Physics,  9,  53  (196?) 

2.  M.  N.  Rosenbluth,  A.  Simon,  Phys .  Fluids  B,    1300  (I965) 
and  9.,  726  (1966) 

3.  G.  Emmert  and  G.  Schmidt,  to  be  published. 
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Normal  Modes  and  Space-Time  Correlations 

J.  K.  FERGUS 

New  York  University 
Courant  Institute  of  Mathematical  Sciences 

G.  J.  YEVICK 

Stevens  Institute  of  Technology 

The  normal  mode  analysis  of  space-time  correlations 
in  a  many-body  system  requires  a  mechanism  for  decay, 
which  may  be  attributed  explicitly  to  velocity  or  current 
density  fluctuations.   The  fashion  in  which  these  may  be 
handled  in  an  ensemble  average  of  the  microscopic  hydro- 
dynamic  equations  is  not  a  priori  transparent.   This 
problem  has  been  investigated  by  a  step-by-step  comparison 
with  the  readily  obtained  solution  of  the  linearized  micro- 
scopic Vlasov  equation.   The  advantage  of  this  framework 
is  that  no  specific  approximation  as  to  the  correlation 
structure  has  to  be  imposed. 
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This  report  was  prepared  as  an  account  of 
Oovernment  sponsored  work.   Neither  the 
United  States,  nor  the  Conimlsslon,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  Implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  Information  contained  In  this  report, 
or  that  the  use  of  any  Information, 
apparatus,  method,  or  proceea  disclosed 
In  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  insulting  from 
the  use  of  any  Information,  apparatus, 
method,  or  process  disclosed  In  this 
report. 

As  used  In  the  above,  "person  acting  on  behalf 
of  the  Commission"  Includes  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  Infor- 
mation pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with 
such  contractor. 
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